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Abstract 

I investigate the possibility of using Coulomb drag to detect a precursor of 
the predicted (but as yet not definitively observed) superfluid transition in 
electron-hole coupled quantum wells. The drag transresisitivity P21 is shown 
to be significantly enhanced above the transition temperature T c and to di- 
verge logarithmically as T — > T c + , due to electron-hole pairing fluctuations 
which are somewhat analogous to the Maki-Thompson contribution to con- 
ductivity in metals above the superconducting T c . The enhancement in P21 is 

estimated to be detectable at temperatures much higher than T c . 
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Under suitable conditions, an attractive interaction between fermions causes a pairing 
instability that leads to phase transition, such as superconductivity, where phonons mediate 
an effective attractive interaction between electrons. Since the direct Coulomb interaction 
is in general stronger than effective phonon-mediated interactions, a coupled electron-hole 
gas might be expected to readily undergo an equivalent phase transition. This possibility 
was theoretically investigated in the sixties, and it was predicted that electron-hole systems 
in bulk systems would form an condensate at low enough temperatures Despite a 

formal similarity between electron-hole condensate and the BCS superconducting states, the 
former was shown to be in fact an insulator. Subsequently, it was proposed || that separated 
electron-hole systems could form a superfluid state, due to the suppresion of transitions that 
fix the phase of the order parameter thereby producing an insulator. The separation of the 
electron and hole components in semiconductor systems also has an important experimental 
advantage; it inhibits the recombination of the electrons and holes which would otherwise 
rapidly deplete the system. 

The advent of fabrication techniques such as molecular beam epitaxy has made it possi- 
ble to produce high quality electron-hole double well systems and to test these predictions. 
Experimental searches for double layer electron-hole condensation have generally used op- 
tical techniques. While there are reports of optical signatures [f|] hinting at a transition, 
there has been no definitive observation of a superfluid condensation. Recently, Vignale and 
MacDonald || proposed the superfluid transition could be definitively identified through a 
completely different route: using Coulomb drag. 

In Coulomb drag experiments, each well in the coupled well system is individually elec- 
trically contacted, making it possible to isolate a driving voltage in just one of the two wells. 
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When an electrical current is driven through one well, the interlayer interactions produce a 
measurable drag signal in the other. The measured quantity is typically the transresistivity 
P21 = E2/J1, where J\ is the current density in the driving layer and E2 is the electric 
field response in the other layer when it in an open circuit. There has been considerable 
experimental and theoretical J7J activity in drag in the past decade. 

For most typical experimental parameters, the interlayer interaction is a weak pertur- 
bation, in which case drag is adequately described within the Born approximation; i.e., the 
carriers in the adjacent wells are uncorrelated, and carriers in one well occasionally col- 
lide with carriers in the other. Typically, p2\ is much smaller than the resistivities of the 
individual layers, pn and p22- Furthermore, in this weak-coupling regime the Born approx- 
imation unambiguously predicts p%\ — > as the temperature T goes to zero, because the 
scattering phase space for the quasiparticles vanishes at zero temperature. What Vignale 
and MacDonald |3j showed is that when electron-hole condensation occurs, p 2 i becomes 
comparable in magnitude to pn and P22, and diverges when T — > 0. As yet, this has not 
been experimentally observed. 

Even when the electrons and holes are not in a condensed state, interlayer correlations can 
affect the transresisitivity f|. In the closely related case of drag of composite fermions (elec- 
trons bound to two magnetic flux quanta), the build-up of interlayer correlations produces 
unique signatures in the transresistivity which differ markedly from the weak-coupling 
result p2i{T —>■ 0) = 0. Deviations from the weak-coupling result in composite fermion 
drag at low temperatures have in fact been reported [RJ . Analogously, deviations from the 
weak-couling form of p2i{T) in electron-hole systems would indicate the presence of inter- 
layer correlations which foreshadow an electron-hole condensation; thus Coulomb drag can 
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be used to "prospect" for situations where electron-hole condensation could occur. 

In this paper, I evaluate the transresistivity in the normal state of an electron-hole double 
quantum well system, going beyond the Born approximation and taking account correlations 
between the two wells. I use the Kubo formalism to calculate 021 = J1/E2, from which the 
transresistivity is given by p 21 ~ — °"2i( "ii°'22) x - The Hamiltonian I use to describe this 
system is 

H = E E^(^)c!(k l )Q(k i ) + 5: / ) 1 (q)p 2 (-q)V(q) + A m p. (1) 
i=l,2 ki q 

Here are the particle field operators of layer i = 1,2, (the drive and drag layers, re- 
spectively), pi is the particle density operator and Hi mp is the contribution from static 



impurities, which are assumed to be uncorrelated "Particle" refers to either elec- 

tron or hole. Within linear response, it does not matter which layer is driven. The static 
transconductivity cr 2 i = — eie 2 lim^o Im[II(u;)/u;] where ti are the charges in layer i and 
Il(if2 n ) = A -1 Jq dr e lUnT (ji(ir)j 2 (0)) is the current-current correlation function (A = 
sample area, Q n = Matsubara boson frequency). In this paper, h,k B = 1. 



To evaluate <r 2 i, the standard diagrammatic techniques |12| , |13| are used. The terms I 
concentrate on are related to the Maki-Thompson terms used to study effects of supercon- 
ducting fluctuations in the conductivity of metals above T c , and the Feynman diagrams are 
shown in Figs. l(a,b). Physically, these terms describe "the effect of ephemeral Cooper pairs 
on the conductivity" [13]]. Composite fermion drag PJTOfl differs from electron-hole drag in 



that the composite fermions have relatively large masses and the phase-breaking processes 
are dominant [[OJ, which favor the Asmalazov-Larkin terms [|TJ],|TIJ. It is assumed in this 
paper that the samples are clean and phase-breaking processes are weak. The correlations 
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are described by T-matrices which are given by the Bethe-Salpeter equation, shown in Figs. 
l(c,d). T pp and T pa correspond to the particle-particle (pp) and and particle-antiparticle 
(pa) channel, respectively (The term "antiparticle" rather than the standard "hole" is used 
to avoid the dual use of the word "hole.") 

If the interlayer interaction is assumed to be static, then the frequency dependence of 
Tpp {Tpa) comes only from the sum (difference) of the energies in the vertices. Then, from 



standard diagrammatic rules [17| , 
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n(zO n ) = -— ^ V lA k l) v 2,x{ k l) F ( k l, k 2; *^n), ( 2a ) 

^ kik 2 



F(ki, k 2 ; iVt n )= T 2 Gifa, ih,n + ifi n )Gi(k 2 , ik ijTl ) G 2(k 2 , ik 2 , n + iQ n )G 2 (k 2 , ik 2>n ) 

{ k pp\Tpp (Ppp] iki tn + ik 2 ^ n + zfi n ) |k pp ) + (kp a |7^, a (P pa ; iki jn ik 2)n ) |k pa ) . (2b) 



Here, Gi is the particle Green function in layer i, and k^ n are Matsubara fermion frequencies. 
The particle velocities are given by v, = kj/mj. For convenience, the momentum depen- 
dences of the T-matrices are parameterized by the center-of-mass coordinates, Ppp = ki±k 2 

pa 

and kw = x 2 ki =p £ik 2 , where Xi = mij{m\ + m 2 ). 

pa 

To obtain the static transresisitivty, the analytic continuation ifl n — > uj + i0 + , followed 
by the limit to — > 0, must be taken. To simplify the expressions, I assume that both wells 
are weakly disordered (justified by the excellent quality of the two-dimensional electron and 
hole gases produced by molecular beam epitaxy), and hence only terms to lowest order in 
Tj~ , the inverse particle lifetime, are kept ||18|| . One obtains after some algebra 

limlm[^(k 1 ,k 2 ^ + z0 + )] ~ £ r ^ r ^ AlfaM) 

n F (uj[) n F (uj' 2 ) n B {-u[ - uj' 2 ){k pp \lm[T pp (P; u[ + lu' 2 -M0 + )|k pp ) 



+n F (u' 1 ) n F (-u)' 2 ) n B (-uj' l + u' 2 ) (k pa |Im[T pa (P; lj[ - lj' 2 + zO + )]|k pa ). 



(3) 



where Ai(k,u) is the spectral function of the particle in layer i and ub and np are the 



bose and fermi functions, respectively. In the weak disorder limit, one can approximate fI9 
A 2 (kj, u) ~ 47TTj(5(a; — £kj> where £ ki is the particle kinetic energy relative to the chemical 
potential in layer i. Substituting this into Eq. (|3|) and using the generalized optical theorem 





(k\lm[Tpp(Ppp,uj + iO + )\k) 



pa pa 



[k\Tpp(P P P,u + zO + )|k + q) 



pa pa 



x 



lim < Im[V G 1 (x 1 P + k, iuo n + iu) n )G 2 (x{P - k, =Rcu n ) > 

^ tu) n ) 



(4) 



yields, from Eqs. (|2a|), (|2b|) and the expression for cr 2 i in terms of II, 



C21 



^7rrir 2 eie2 /" <iki /" c?k 2 



T 



(2vr) 2 y (2n) 2 J (2tt) 2 

nir(-^ k2 _q) |(kpp|Tpp(Ppp, f kj + £k 2 )|kpp + q)| 5(Ck! + fk a - fkx+q - £k 2 -q) - (fk x )«F (~& a ) 
^F(-^ki+q) ^F(^k 2 +q) |(kp a |Tp a (Pp a ,^ kl - £ k2 )|kp a + q)| 2 5^ki - £k 2 ~ {£ki+q ~ £k 2 +q}) j- (5) 

Eq. (jjD can to a certain extent be interpreted within the framework of the semiclassical 
Boltzmann equation |2]J . The first terms within the integral of Eq. (^) corresponds to the 
contribution of particle-particle scattering from state k 1; k 2 with exchange of momentum 
q, when the distribution of layer 1 is shifted by the driving electric field. The second term 
corresponds to the contribution due to Pauli blocking of electrons going into state k x + q, 
or alternatively, one can look at it as scattering of antiparticles in the Fermi liquid in layer 
1 from ki + q with particles in layer 2. Note that with a purely semiclassical argument, it 
is impossible to have different scattering rates T pp for particle-particle and T pa for particle- 
antiparticle collisions, as is obtained above from a full many-body calculation. In the Born 
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approximation limit, (k|T pp (P, u> + i0 + ) |k + q) = (k|T pa (P, u + iO + )\k + q) = V(q,u), the 



screened interlayer interaction, one regains the standard weak-coupling result [P,|I5,pT-25 
which typically involves a product of the individual layer susceptibilities (or, in the general 
case, a closely related function jnfl). The above shows that when correlations between layers 
are included, the expression for o" 2 i cannot be written in the standard weak-coupling form. 

When the T-matrix is regular, temperature dependence of the transresistivity is not 
qualitatively different from the weak-coupling case. For example, for hard-core fermions, 
as T — > the T-matrix remains finite and consequently the transresistivity vanishes, in 
spite of the presence of interlayer correlations. This result is consistent with the claim that 
a nonzero transresistivity at T = is possible only with a phase transition [23|. When a 



phase transition such as superconductivity or electron-hole condensation takes place, the 
T-matrix acquires a singularity at the (mean-field) transition temperature T c [25[], below 
which the T-matrix approximation is invalid. Below, I examine 021 for T > T c . 

For simplicity, to obtain qualitative features let us assume a simple local interlayer in- 
teraction, V(q) = Vq < 0. Then the T-matrices within the Bethe-Salpeter approximation 
are dependent only on the total momentum P and energy uj of the incoming particles. The 
T pp (P,c<j) channel diverges at T c , u = and P = 0. Its explicit form is 

T PP {P^-T) = - ^\ (6a) 
Xpp{P,v;T) = / -— - — — ■ 6b 

Here, E k = k 2 {ra{ 1 +m 2 l )/2 and E P = P 2 /2(m 1 +m 2 ). 

The singularity in the T pp occurs when \ pp = — |Vo| -1 . This occurs at the highest tem- 
perature when P = 0, u = and kp,i = &f,2 = kp (matched fermi surfaces). Expanding 



Xpp(T) about P = 0, ijj = and $ = [k 2 F2 ~ kp 1 ]/2(mi + m 2 ) gives 

T PP (P, u; T) = (7) 

where p red = m 1 m 2 /2Ti{m 1 + m 2 ), <5x(T) = 21og(T/T c ), a w = 7r/(4k B T), a$ = 
7C(3)/(2T 2 tt 2 ) « 0.43/T 2 and a P = 7((3)£ kF /(T 2 n 2 ) « 0.85£ fcF /T 2 |2§. Substituting 
Eq. (0) into Eq. (H) yields (assuming a 21 <C crn<722, and reintroducing h) 

h T 2 87T 2 / ; \ 
P21 = -^^7C(3T J ( Xl ' $ ' S ) (8a) 

j^, *, s) = r * tj^im _ (8b) 

J-oo y [cosh(y/2) + cosh(y(l/2 - Xl ) + $)] 2 
where s = 811- 1 log (£) + 14C(3)tt- 3 <I 2 and <£ = $/T c . 

For matched electron and hole densities (i.e., $ = 0), 

( l/s if s > 1, 

X(x lt 0, s = 87T" 1 log(T/T c )) ~ (9) 

1 mm ifs<l, 
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ignoring prefactors which depend weakly on 2^ and are of order 1. Thus, the transresistivity 
diverges logarithmically as T — > T c + . Furthermore, p 2 i(T 3> T c ) ~ l/log(T/T c ), and this 
relatively slow fall-off with increasing temperature implies that the electron-hole pairing 
fluctuation enhancement of p 2 i can be seen well above T c . A comparison of Eq. (jBaD with the 
weak-coupling result |2]J P2i,weak — — hT 2 ((3)iimim 2 / (8h 2 eie 2 k F1 k F2 q TF: iq TF!2 d 4 ), implies 



that the pairing fluctuation contribution to the transresistivity is bigger than p 2 i iWea k when 

log ^— J < 7 ^(3)]2 a:i ( 1 ~~ x i) k F<lTF,iqTF,2d 4 . (10) 

where gxF,i = 2rriie 2 /Kh 2 (k = dielectric constant) is the Thomas- Fermi screening length 
in layer i and d is the center-to-center well separation. For GaAs parameters, with kp = 
10 6 cm -1 and d = 4x 10~ 6 cm, the number on the right of Eq. (0) is approximately 5 x 10 3 . 
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Clearly, this result should not be interpreted quantitatively, since several approximations 
and assumptions have been used in this calculation. The local interlayer interaction no 
doubt leads to a significant overestimate of the pairing fluctuation contribution, since it fails 
to cut off the large momentum transfer contributions, unlike the more realistic interlayer 
Coulomb interaction. Furthermore, the above calculation does not take into account other 
factors that tend to impede electron-hole condensation. These include uncorrelated impurity 
potentials in the electron and hole layers (negating the the time-reversal argument which 
makes non-magnetic impurity scattering irrelevant in superconductivity) and band-structure 



effects |[27|| . Nevertheless, the above result strongly suggests that the enhancement can be 
experimentally detected. 

For umatched electron-hole densities (i.e., $ 7^ 0), the density dependence of the pairing 
fluctuation enhanced P21 differs from the weak-coupling result. For kp,i ~ kp^, the weak- 
coupling transresitivity goes as P2i,weak ~ (^f,i^f,2) -3 ; i.e. there is a monotonic increase with 
decrease of either density. On the other hand, since the nesting of the electron and hole fermi 
surfaces is a necessary condition for electron-hole condensation, in the pairing fluctuation 
enhanced case the transresistivity peaks at matched densities, p2i(& = 0) — P2i( < £) oc (<1>/T) 2 . 
It should be noted that a peak in the P21 at Izf,i = &f,2 is also the signature of a phonon 
mediated interaction. However, the phonon mediated interaction falls rapidly (~ T 6 ) below 



the Bloch-Griineisen temperature |28| which is typically a few degrees K for GaAs, so if a 
the peak is observed at low enough temperatures the phonon mechanism can be ruled out. 

To conclude, I have presented a strong-coupling theory of Coulomb drag, within the T- 
matrix approximation. Applying this theory to electron-hole double quantum well systems, 
I find that the pairing fluctuations lead to a large enhancement in the transresistivity above 



the electron-hole condensation temperature T c , and this effect could be used to identify 
promising candidates for the observation of electron-hole condensation. 

This work was initiated at Mikroelektronik Centret, Danmarks Tekniske Universitet, 
and partially supported by a University of Akron Summer Research Fellowship. I gratefully 
acknowledge useful discussions with Martin Chr. B0nsager, Karsten Flensberg, Antti-Pekka 
Jauho and John W. Wilkins. 
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FIG. 1. Feynman diagrams used in calculating the transconductivity. Figs, (a) and (b) show 
the particle-particle and particle-antiparticle channel contributions to 021 • The black dots are the 
current vertices, the numbers above them denote the layer indices, and the arrows are the particle 
Green functions. The T pp and T pa are the particle-particle and particle-antiparticle T-matrices 
which are calculated using the Bethe-Salpeter equation, shown diagramatically in Figs, (c) and 
(d), respectively. 
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